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ABSTRACT 

We present a detailed analysis of a large-scale galactic outflow in the CGM of a massive (Mh ~ Mq), 

star forming (~ 6.9 Mq yr”'^), sub-E, (~ 0.5Lg) galaxy at 2 ; = 0.39853 that exhibits a wealth of metal-line 
absorption in the spectra of the background quasar Q 0122 — 003 at an impact parameter of 163 kpc. The 
galaxy inclination angle {i = 63°) and the azimuthal angle (fl) = 73°) imply that the QSO sightline is passing 
through the projected minor-axis of the galaxy. The absorption system shows a multiphase, multicomponent 
structure with ultra-strong, wide velocity spread O VI (logTV = 15.16 ± 0.04, Augo = 419 km s”'^) and N V 
(log N = 14.69 ± 0.07, Augo = 285 km s”'^) lines that are extremely rare in the literature. The highly ionized 
absorption components are well explained as arising in a low density (~ 10“^ ^ cm“^), diffuse (^ 10 kpc), 
cool (^ 10^ K) photoionized gas with a super-solar metallicity ([X/H] > 0.3). Erom the observed narrowness 
of the Ly/3 profile, the non-detection of S IV absorption, and the presence of strong C IV absorption in the 
low-resolution EOS spectrum we rule out equilibrium/non-equilibrium collisional ionization models. The low- 
ionization photoionized gas with a density of ^ 10“^'^ cm“'^ and a metallicity of [X/H] > —1.4 is possibly 
tracing recycled halo gas. We estimate an outflow mass of ~ 2 x lO'^® Mq, a mass-flow rate of ~ 54 Mq yT~^, 
a kinetic luminosity of ~ 9 x 10''^ erg s”'^, and a mass loading factor of ^ 8 for the outflowing high-ionization 
gas. These are consistent with the properties of “down-the-barrel” outflows from infrared-luminous starbursts 
as studied by Rupke et al. Such powerful, large-scale, metal-rich outflows are the primary means of sufficient 
mechanical and chemical feedback as invoked in theoretical models of galaxy formation and evolution. 

Subject headings: galaxies:formation - galaxiesihaloes - quasarsiabsorption lines - quasarundividual 
(Q 0122-003) 


1. INTRODUCTION 

Accretion of gas from the intergalactic medium (IGM) to 
the interstellar medium (ISM) and large-scale galactic out¬ 
flows are the two principal components of current theoret¬ 
ical models of galaxy formation and evolution (Springel & 
Hernquist 2003; Keres et al. 2005; Dekel et al. 2009; Oppen- 
heimer et al. 2010; Dave et al. 201 la,b, 2012). Gas accretion 
is thought to be bimodal: (a) the “hot mode” in which infalling 
gas onto a massive halo (Mh > IO'^^Mq) gets shock-heated 
near the halo virial radius (i?vir) to the halo virial temperature 
(Tvir ^ 10® K) and (b) the “cold mode” in which gas enters 
well inside i?vir in a less massive halo with T ^ 10^— 10® K, 
without any subsequent shock heating (Keres et al. 2005; 
Dekel & Birnboim 2006; Dekel et al. 2009; Keres et al. 2009). 
Large-scale wind/outflows, on the other hand, are the primary 
mechanism by which energy, baryons, and metals are recy¬ 
cled in galaxies and transported to the IGM (see Veilleux et al. 
2005, for a review). Such mechanical and chemical feedback 
is essential to regulate star-formation rates (SERs) and stel¬ 
lar masses (M*) in theoretical models (Springel & Hernquist 
2003; Dave et al. 2011b). However, it has been very chal¬ 
lenging to directly detect these processes and thereby test the 
theoretical models. 

The infall of metal-poor gas and the outflow of metal- 
enriched gas both take place through the circumgalactic 
medium (CGM). Understanding the physical conditions and 
chemical abundances of the CGM, thus, is of utmost impor¬ 


tance in the context of how galaxies form and evolve with 
cosmic time. In particular, the amount of diffuse baryons 
and metals, their geometrical distributions, thermal/ionization 
states, and the kinematics of the gas in the CGM provide im¬ 
portant clues on how galaxies acquire their gas and how they 
recycle it in the form of feedback. Ultraviolet (UV) absorp¬ 
tion lines originating from the CGM of a foreground galaxy 
in the spectrum of a nearby background quasar (QSO) act as 
sensitive probes of the physical and chemical conditions of 
the CGM gas. A significant amount of work has been done, in 
the recent past, to better understand the CGM gas (e.g. Chen 
et al. 2001; Chen & Mulchaey 2009; Prochaska et al. 2011; 
Kacprzak et al. 2008, 2010, 2012; Tumlinson et al. 2011b; 
Thom et al. 2012; Nielsen et al. 2013a,b; Werk et al. 2013; 
Stocke et al. 2013; Mathes et al. 2014; Werk et al. 2014; Bor- 
doloi et al. 2014a; Borthakuret al. 2015). 

The primary goals of these above studies are to under¬ 
stand the covering fractions of different phases of the CGM 
gas and their distributions around different galaxy types, i.e., 
star-forming (SE) and non-SE. Eor example, in the recent 
COS-Halos survey it has been seen that the covering frac¬ 
tion of Ovi absorption, with rest-frame equivalent width 
Wr > 0.1 A, is significantly higher around SE galaxies 
as compared to non-SE galaxies (Tumlinson et al. 2011b). 
The covering fractions of low-ions (e.g. Mg II, Si III) and 
H I, on the other hand, do not show any preference on the 
galaxy types (see for example Eigure 12 of Werk et al. 2013). 
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Nonetheless, Mg II absorbers are known to exhibit a strong 
dependence on galaxy orientations (Bordoloi et al. 2011; 
Kacprzaket al. 2011,2012; Bouche et al. 2012; Bordoloi et al. 
2014b). For example, Kacprzak et al. (2012) have reported 
a bimodality in the azimuthal angle distribution of halo gas 
traced by Mg II absorption, such that it prefers to exist near 
the projected major and minor axes around SF galaxies. Non- 
SF galaxies, on the contrary, do not show any such preference. 
The kinematics of Mg II absorbers are further found to depend 
on galaxy orientation. For example, in a recent work, Nielsen 
et al. (2015) have shown that face-on star forming galaxies 
viewed along their projected minor-axis have the largest ve¬ 
locity dispersion, suggesting that some fraction of Mg II ab¬ 
sorbers directly traces bi-conical outflows. Therefore, con¬ 
necting absorption properties with both star formation rate 
of host-galaxy and its orientation relative to the QSO sight¬ 
line are essential in order to comprehend the complete picture 
of the mechanisms by which baryons are processed through 
galaxies, i.e., the baron cycle. 

One of the important findings of the COS-FIalos survey is 
that O VI is omnipresent in the halos of SF galaxies which har¬ 
bor a major reservoir of galactic metals (i.e. Tumlinson et al. 
201 lb). Outflows from the SF galaxies are thought to be the 
origin of highly ionized oxygen in the COM. However, the 
outflows need not be active at the time of observations. Active 
outflows that are detected primarily via low-ions (e.g. Na I, 
Fell, Mg II) absorption in the spectra of host-galaxies (i.e., 
the “down-the-barrel” outflows) are also found to be ubiq¬ 
uitous at both high and low redshifts (Shapley et al. 2003; 
Veilleux et al. 2005; Martin 2005; Rupke et al. 2005; Rubin 
et al. 2010,2014). Despite the ubiquity of winds, several basic 
questions that are key to quantifying galaxy feedback remain 
unanswered; How far do they propagate? What is the baryon 
mass that are processed through them? What is the mass- 
flow rate through winds? What are the mass loading factors 
and kinetic power of the winds? To what degree are metals 
processed through them? This is essentially because the loca¬ 
tion of the outflow with respect to host-galaxy is an unknown 
in “down-the-barrel” outflows. Galaxy outflow probed by a 
background QSO has the advantage of having the minimum 
distance between host-galaxy and outflowing material. For 
example, by analyzing a post-starburst outflow from a galaxy 
at z ~ 0.927 detected in the spectra of QSO PG 1206H-459 
at an impact parameter of ~ 68 kpc, Tripp et al. (201 1) have 
shown that the entrained gas mass could be as large as ^ few 
100 X 10®M 0 . Numerous strong and large velocity spread 
(Au > 400 km s“^) metal absorption lines along with so¬ 
lar to super-solar metallicities in different absorption compo¬ 
nents suggest that the absorber is tracing an active outflow. 
Such wide velocity spread O VI absorbers at low-z are also 
reported by Tumlinson et al. (2011a) and Muzahid (2014). In 
both these cases metallicity of the O VI bearing gas is low (e.g. 
[X/H] ^ —1.0) which suggest that Ovi is probably tracing 
an “ancient outflow” (Ford et al. 2014) rather than an active 
wind. Interestingly, strong N V is detected in PG 1206H-459 
which is not present in the latter two systems. We note that the 
location of the absorbing gas with respect to the host-galaxy’s 
projected major and minor axes are not known for any of these 
strong O VI absorbers. 

Here we present a detailed study of an active outflow in 
the CGM of a star-forming (SFR ^ 6.9Mq yr“^ ), sub-L* 
(- 0.5L|j), Sbc galaxy at redshift z = 0.39853. The CGM of 
the galaxy produces several low- and high-ionization absorp¬ 
tion lines in the FLT(UVES) and HST(COS, FOS) spectra 


of the background quasar Q 0122 — 003 at an impact param¬ 
eter of 163 kpc. The observed inclination angle, i = 63° 
and azimuthal angle, 4> = 73° indicates that the QSO sight¬ 
line is at ideal location to probe a large-scale bipolar outflow 
from the host-galaxy (Kacprzak et al. 2012; Bordoloi et al. 
2014b; Fox et al. 2015). This article is organized as follows; 
In Section 2 we present observations of the background QSO 
and the host-galaxy. Analysis of the absorption system, the 
host-galaxy, and overall photoionization (PI) models, based 
on the observed total column densities, are presented in Sec¬ 
tion 3. Detailed component-by-component ionization models 
are presented in Appendix B for the interested readers. In 
Section 4 we discuss the implications of our findings. The 
results of our study are summarized in Section 5. Through¬ 
out this article, we adopt a flat ACDM cosmology with Hq = 
71 km s“^Mpc“^, Hm = 0.3, and Ha = 0.7. Abundances are 
given in the notation [X/Y] = log(X/Y) — log(X/Y )0 with 
solar abundances taken from Asplund et al. (2009). Atomic 
data for various elements are taken from Morton (2003). All 
the distances given are proper (physical) distances. 

2 . OBSERVATIONS 

2.1. Absorption Data 

2.1.1. HST/COS 

UV spectra of the background QSO Q 0122-003 were ob¬ 
tained using HST/COS Cycle-21 observations, under pro¬ 
gram ID; GO-13398 as a part of our “Multiphase Galaxy Ha¬ 
los” survey. The properties of COS and its in-flight opera¬ 
tions are discussed by Osterman et al. (2011) and Green et al. 
(2012). The observations consist of G160M far-UV (EUV) 
grating integrations with a wavelength coverage of 1410- 
1780 A at a medium resolution of i? ~ 20,000 (EWHM ^18 
km s“^). The data were retrieved from the HST archive and 
reduced using the STScI CALCOS v2.21 pipeline software. 
The reduced data were flux calibrated. To increase the spec¬ 
tral S/N, individual G160M integrations were aligned and 
coadded using the IDE code “coadd^xld” developed by Dan- 
forth et al. (2010)'. The combined spectrum has a S/N ~ 5 - 
10 per resolution element. As the COS EUV spectra are sig¬ 
nificantly oversampled (i.e. six raw pixels per resolution ele¬ 
ment), we binned the data by three pixels. This improves S/N 
per pixel by a factor of v/3- All our measurements and anal¬ 
yses were performed on the binned data. Continuum normal¬ 
ization was done by fitting the line-free regions with smooth 
low-order polynomials. 

2.1.2. ULT/UVES 

The optical spectrum of Q 0122-003 was obtained with the 
Ultraviolet and Visible Echelle Spectrograph (UVES; Dekker 
et al. 2000) mounted on the European Southern Observatory 
Kueyen 8.2-m telescope at the Paranal Observatory during 
28-31 July, 2005 under program ID; 075.A-0841. The UVES 
observations covered 3290-9466 A at spectral resolution of 
R ~ 45,000 (EHWM ~ 6.6 km s“'; 3 pixels per resolution 
element). The spectrum has the highest S/N of ~ 20 per 
pixel roughly between 5000-7000 A. The S/N drops consid¬ 
erably to ~ 5 per pixel in both the blue and red ends of the 
spectrum. We refer the reader to Kacprzak et al. (2011) and 
Evans (2011) for the full details of data reduction procedure. 

* http://casa.colorado.edu/danforth/science/cos/costools.html 
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2.2. Galaxy Data 

A 2100 second HST/WFPC2 F702W image of the 
Q 0122-003 field was obtained under program ID; 6619. The 
reduced and calibrated image was obtained from the WFPC- 
2 Associations Science Products Pipeline (WASPP). The im¬ 
age has a limiting magnitude of 26, which translates to an 
Mb = -14.5 and L = 0.002L* at z = 0.4. 

Spectra of seven galaxies in the Q 0122-003 field were ob¬ 
tained using the Keck Echelle Spectrograph and Imager (ESI; 
Sheinis et al. 2002) on 2014 December 13. The mean seeing 
was 0.8" (FWHM) with clear skies. Exposure times range 
between 1800-3600 seconds. The slit is 20" long and 1" wide 
and we used 2x2 on-chip CCD binning. Binning by two in 
the spatial directions results in pixel sizes of 0.27 — 0.34" 
over the echelle orders of interest. The wavelength cover¬ 
age of ESI is ~ 4000-10,000 A, which provides coverages 
of multiple emission lines such as Oil doublets, H/3, Olll 
doublets, Ha, [N ll] doublets with a velocity dispersion of 
22 km s“^ pixel“^ when binning by two in the spectral di¬ 
rection (FWHM ~ 90 km s“^). The spectra were reduced 
using the standard echelle package in IRAF along with stan¬ 
dard calibrations and were vacuum and heliocentric velocity 
corrected. Spectra were flux calibrated using the standard star 
Feige34 taken during the night of the observation. We have 
made no corrections for slit loss or Galactic reddening. 

3. ANALYSIS 
3.1. Absorption Analysis 

The Lya, Ly/3 and different ionic transitions originating 
from the Zahs = 0.39853 absorber are shown in Figure 1 . Mg II 
absorption, covered by the high resolution UVES spectrum, 
is detected in three velocity components. We refer to them 
as L-1 (v ^ —10 km s“^), L-2 (v ^ -1-180 km s“^), and 
L-3 (v ^ -1-200 km s“^). Mg I absorption is detected only 
in the two strong Mg II components at u ^ -1-200 km s“^ 
(i.e. in L-2 and L-3). A weak Ee II absorption is also de¬ 
tected in these two components. Using the VPFIT^ software 
we fit Mg I, Mg II, and Ee II lines simultaneously assuming the 
6 -parameters of magnesium and iron are connected via pure 
thermal broadening. The Voigt profile fit parameters are sum¬ 
marized in Table 1. In the cases of non-detections, we present 
standard 3cr upper limits estimated from the error spectrum. 

As we use only the G160M grating data, any transi¬ 
tions with rest-frame wavelength below 1010 A and above 
1270 A are not covered by the COS spectrum. The detected 
C II A1036 and N II A1083 absorption clearly trace the Mg II 
line kinematics seen in the UVES spectrum. However, the 
component structure at u ~ -1-200 km s“^ is not resolved 
due to the lower spectral resolution of COS. Note that the 
N II A1083 absorption at u ^ 0 km s“^ is blended with 
Si III A1206 absorption from another known intervening sys¬ 
tem at Zabs = 0.25647. The Sill A1193 is blended with the 
Galactic Al II A1670 absorption. The Si II A1190, A1260 and 
Si III A1206 transitions are too noisy to get reasonable fitting 
parameters. Moreover, the latter two transitions are partially 
blended by unrelated absorption. As the component structure 
is not resolved and the data are noisy, except for the weak 
C II component at ^ 0 km s“^, we do not attempt to fit the 
low-ions detected in the COS spectrum. Instead, we gener¬ 
ate synthetic profiles corresponding to our adopted PI model 
(Appendix B.2) and compare with the observed data. 

^ http://www.ast.cam.ac.uk/ rfc/vpfit.html 


We note here that the line spread function (LSE) of the 
COS spectrograph is not a Gaussian. A characterization of the 
non-Gaussian LSE for COS can be found in Ghavamian et al. 
(2009) and Kriss (2011). Eor our Voigt profile fit analysis we 
adopt the latest LSE given by Kriss (2011). The LSE was ob¬ 
tained by interpolating the LSE tables at the observed central 
wavelength for each absorption line and was convolved with 
the model Voigt profile while fitting absorption lines or gen¬ 
erating synthetic profiles using VPFIT. Since there is no Lya 
forest crowding at low-z and the interested absorption lines do 
not fall on top of any emission line, we do not account for the 
continuum fitting error in the absorption line fit parameters. 

The Ovi AA1031,1037 absorption is detected over a ve¬ 
locity spread of > 500 km s“^. A formal measurement gives 
Augo = 419 km s“^ with an optical depth-weighted absorp¬ 
tion redshift of z = 0.39909 (see e.g. Eigure 3 of Muzahid 
et al. 2012). As discussed in Section 4, this is among the four 
largest velocity spread intervening O VI absorbers at z < 1 
known to date. A minimum of five components are required 
to fit the Ovi doublets satisfactorily. We refer to them as 
H-1 through H-5 from lowest velocities to highest veloci¬ 
ties. Note that the O VI A1037 absorption in the H-1 and H-2 
components is blended with the C II A1036 absorption. The 
apparent column density profiles do not suggest any hidden 
saturation in the three redward O VI components (i.e. H-3, 
H^, and H-5) which are not affected by the C II contami¬ 
nation. It is intriguing to note that, unlike most low-z O VI 
absorbers, strong N V is detected in these three components. 
However, the N V profiles appear noisy and the N V A1242 
line is partially blended. To obtain reasonable fitting param¬ 
eters we fit both O VI and N V simultaneously assuming pure 
thermal broadening. A pure non-thermal assumption also pro¬ 
vides similar component column densities. Voigt profile fit 
parameters for these components are given in Table 2. In the 
cases of non-detections we provide standard 3(7 upper limits. 

It is apparent from the Eigure 1 that the absorption kine¬ 
matics of Mg II and O VI are significantly different. Mg II 
shows two distinct absorption clumps at u ^ —10 km s“^ 
and at u ^ -1-200 km s“^ with Augg of ^ 25 and ~ 28 
km s“^, respectively. Ovi, on the contrary, shows a con¬ 
tiguous absorption with a Augg of 419 km s“^. No low- 
ions are detected in the H-1, H-3, and H-5 components. 
Mg II is detected at velocities similar to those of the H-2 and 
H-4 components. However, the O VI components show much 
larger 6-parameters as compared to the Mg II components. 
The data, thus, clearly indicate that the absorbing gas has mul¬ 
tiple phases with different densities (and/or different temper¬ 
atures) giving rise to different absorption kinematics. 

Here, we refer the reader to Appendix A for a detailed de¬ 
scription of how we estimate A(H l) in the different high- and 
low-ionization absorption components. In brief, we estimate 
log A(Hl) in the range 18.3-18.6 and 14.6-15.3 for the low- 
and high-ionization components, respectively. 

3.2. Galaxy Analysis 

The HST/WFPC2 E702W image of the Q 0122-003 field 
is shown in the left panel of Eigure 2. There are seven galax¬ 
ies with confirmed spectroscopic redshifts (see Table 3 for de¬ 
tails). Galaxy apparent Vega-magnitudes were determined us¬ 
ing L5cr isophotes from Source Extractor (Bertin & Arnouts 
1996). The sky orientation parameters, such as inclination an¬ 
gles (/) and azimuthal angles ($), for all these galaxies were 
measured using GIM2D (Simard et al. 2002) following the 
methods of Kacprzak et al. (2011). Here we define the az- 
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Figure 1. Velocity plots of different transitions covered in the HST /COS and VLT/UVES spectra. The zero velocity corresponds to the redshift of the host- 
galaxy, i.e. ZgoX = 0.39853. The best fitting Voigt profiles for high- (low-) ions are shown in smooth red (blue) curves on top of data. In the Lya and Ly/5 panels 
the blue and red curves indicate the contributions of low- and high-ionization gas phases, respectively (Appendix A). The magenta cuiwes represent synthetic 
profiles corresponding to our adopted PI models for the low-ions (Appendix B.2), since the COS spectrum is noisy and does not have enough resolution to resolve 
the nari'ow lines. Error in the normalized fiux is shown in green histogram. The vertical blue dotted lines represent the line centroids of the Mg II components. 
The small red ticks in each panel show the line centroids of the O VI components. Unrelated absorption are marked by ‘x’ symbols. Note that O VI A1037 model 
subtracted spectnam is shown in the C II A1036 panel. 


imuthal angle such that at 4) = 0° the quasar line-of-sight lies 
along the galaxy projected major axis, and at 4> = 90° it lies 
along the projected minor axis. 

We use our own fitting program (FITTER: see Churchill 
et al. 2000), which computes best fit Gaussian amplitudes, 
line centers, and widths, in order to obtain emission-line red- 
shifts and line fluxes. The galaxy redshifts as listed in Table 3 
have accuracy ranges from 6-15 km s“^. Here we focus 
our study on the z^ai = 0.39853 galaxy and will present de¬ 
tailed analysis on the other galaxies in future works. We show 
them here to demonstrate the spectroscopic completeness of 
the field. Note that no other galaxy redshift matches the ab¬ 
sorption redshift within Av > ±1000 km s“^, suggesting 
that the host-galaxy is isolated. The host-galaxy has an im¬ 
pact parameter D = 163 kpc, i = 63° and 4> = 73°. 


The SFR of the host-galaxy is estimated from the Ha lu¬ 
minosity using the relation of Hao et al. (2011). The mea¬ 
sured Ha flux of 2.69 x 10“^^ erg cm“^ s“^, corresponding 
to a luminosity of 1.50 x 10^^ erg s“^, leads to a SFR of 
6.9 Mq yr“^. Using the effective radius (3.5 kpc) and el- 
lipticity (0.56) obtained from the GIM2D model we derive 
a star-formation rate density of S* = 0.4 Mq yr”^ kpc“^. 
This is a factor of 4 higher than the known E* threshold 
for driving galactic-scale outflows in local starbursts (Heck¬ 
man 2003). The observed Olll/H/? and Nll/Ha emission 
line ratios (i.e. —0.323 and —0.391 dex respectively) sug¬ 
gest that while most of the emission is due to star formation, 
some AGN contamination could also be present (Baldwin 
et al. 1981). Thus, strictly speaking, the measured SFR and 
E* should be considered as upper limit. We compute a gas- 
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Table 1 

Low-ionization metal line fit parameters. 


Ion 

■^abs 

ID 

6 (km s ^) 

log Af (cm“^) 

Me II 

0.398479±0.000003 

L^l 

9.1 ±0.8 

12.44±0.03 

Cii 



9.1 

13.90±0.14 

Fell 



9.1 

<12.5 

Call 



9.1 

<11.3 

Mgl 



9.1 

<11.2 

Nal 



9.1 

<11.6 

Hi 



9.1“ 

<18.6 

Mg II 

0.399384±0.000002 

L-2 

4.8T0.5 

13.24±0.08 

Fell 



3.1T0.5 

12.66±0.10 

Call 



4.8 

<11.2 

Mgl 



4.8T0.5 

11.62±0.05 

Nal 



4.8 

<11.6 

Hi 



4.8“ 

<18.3 

Mg II 

0.399449±0.000003 

L-3 

6.7±0.6 

13.11±0.04 

Fell 



4.4±0.6 

12.56±0.10 

Call 



6.7 

<11.3 

Mgl 



6.7±0.6 

11.53±0.06 

Nal 



6.7 

<11.4 

Hi 



6.7“ 

<18.3 


Notes - All column density upper limits are quoted at 3cr as estimated from 
the eiTor spectrum for the adopted b-parameter. “The 6-values are too low 
for the gas temperature we derived under PI equilibrium conditions (i.e. T ~ 
10^ K). However, we found that the Af(Hl) limits remain unchanged for 6- 
parameters of 14-20 km s“^ as suggested by our PI models in Appendix B. 

phase oxygen abundance for the host-galaxy using the N2 re¬ 
lation of Pettini (2004) where 12H-log(O/H)=8.90H-0.57xN2 
(N2 = log(Nll/Ha)). 

The galaxy halo mass M h was obtained by halo abundance 
matching. We used the methods described in Churchill et al. 
(2013a). Systematics and uncertainties in the method are 
elaborated in Churchill et al. (2013b). The halo mass library 
is drawn from the Bolshoi A^-body cosmological simulation 
(Klypin et al. 2011) and is matched to the observed r-band 
luminosity function from the COMBO-17 survey (Wolf et al. 
2003). The galaxy r-band and i?-band absolute magnitudes, 
Mr and Mb (respectively), were determined by iT-correcting 
(e.g., Kim et al. 1996) the HST/WFPC2 F702W observed 
magnitude using the Coleman et al. (1980) spectral energy 
distributions (SED) from Bolzonella et al. (2000). The K- 
correction is fully described in Nielsen et al. (2013b). The 
galaxy B-band luminosity, i.e. Lb/L*^ = 0.5, was calcu¬ 
lated by using the linear fit to Mg with redshift from Faber 
et al. (2007). We have no color information for the galaxy. 
Given the late-type appearance of the galaxy morphology in 
the WFPC2 image (see the top-right panel of Figure 2), we 
adopt the Sbc SED for the iC-correction. The adopted halo 
mass is logMh/M© = 12.48 ± 0.16. The iT-correction 
varies by no more than 0.46 from an Irr to an E SED, which 
translates to a halo mass range of log Mji/Mq = 12.4-12.7, 
respectively, across the range of “normal” galaxy SEDs. 

The galaxy virial radius, i?vir, is obtained using the formal¬ 
ism of Bryan & Norman (1998). We obtained i?vir = 278lp 
kpc, where the uncertainty accounts for the uncertainty in 
the adopted Mh. For the line-of-sight impact parameter, we 
probe the CGM of this galaxy with respect to the virial radius 
at the projected location of D/Rvir = 0.6. 

3.3. Photoionization modeling 

In order to understand the physical conditions and the 
chemical abundances in the absorbing gas we run grids of 


Table 2 

High-ionization metal line fit parameters. 


Ion 

-^abs 

ID 

6(kms-^) 

log N (cm ) 

Ovi 

0.398153± 0.000023 

H-1 

28.3±6.3 

14.26± 0.09 

N V 



28.3 

<13.7 

Siv 



28.3 

<14.0 

Hi 



28.3 

<14.9(15.04±0.09) 

Ovi 

0.398481± 0.000026 

H-2 

35.4±10.5 

14.28± 0.09 

N V 



35.4 

<13.7 

Siv 



35.4 

<14.0 

Hi 



35.4 

<15.2(15.36±0.11) 

Ovi 

0.398913± 0.000009 

H-3 

32.5± 3.2 

14.62± 0.03 

N V 



34.7± 0.0“ 

14.22± 0.08 

Siv 



32.5 

<14.1 

Sim 



32.5 

<12.8 

Hi 



32.5 

<14.8(14.83±0.05) 

Ovi 

0.399443± 0.000009 

H=i 

42.4± 3.6 

14.61± 0.02 

N V 



45.3± 0.0“ 

14.28± 0.06 

Siv 



42.4 

<14.2 

Hi 



42.4 

<15.3 (15.35±0.05) 

Ovi 

0.399964± 0.000010 

H-5 

32.5± 3.2 

14.38± 0.03 

N V 



34.8± 0.0“ 

14.12± 0.08 

Siv 



32.5 

<14.1 

Sim 



32.5 

<12.8 

Hi 



32.5 

<14.6(14.64±0.04) 

Civ'’ 

(AOD measurements) 


191 

>14.7 


Notes - All column density upper limits are quoted at 3fT as estimated from 
the error spectrum for the adopted 6-parameter. The Al(Hl) values in the 
parenthesis are obtained assuming that the Ly/3 profile is entirely due to the 
high-ionization absorption components (see Appendix A) with 6-parameters 
tied with 6(0 Vl). “6(N v) are tied with 6(0 Vl) via thermal broadening. 

Apparent Optical Depth (AOD) measurements using FOS spectrum. 

PI models using CLOUDY (vl3.03, last described by Ferland 
et al. 2013). In these models absorbing gas is assumed to 
be a plane parallel slab, exposed to the extra-galactic UV 
background radiation at z = 0.39 (as computed by Haardt 
& Madau 2012) from one side. Here we do not consider 
the effect of a galaxy/stellar radiation field, as it is negligi¬ 
ble at this redshift at a large separation (> 100 kpc) from 
bright (^ L*) galaxies (e.g. Narayanan et al. 2010; Werk et al. 
2014). We note that the host-galaxy for the present system has 
L = 0.5Lg with an impact parameter of 163 kpc. In our mod¬ 
els, the relative abundances of heavy elements are assumed to 
be solar as in Asplund et al. (2009). 

Note that here we estimate the average ionization condi¬ 
tions and abundances for the high- and low-ionization gas us¬ 
ing the total column densities of relevant ions and H I. To¬ 
tal column density is derived by summing the component 
column densities presented in Table 1 and 2. We refer the 
interested readers to Appendix B for a detailed component- 
by-component ionization modeling for both the high- and 
low-ionization phases. In Appendix B we have explored 
both PI and collisional ionization equilibrium (CIE) and non¬ 
equilibrium (non-CIE) models for the high-ionization phase. 
Moreover, we ruled out both the CIE and non-CIE scenarios 
for the O VI bearing high-ionization gas using several obser¬ 
vational constraints such as; (i) the Ly/3 profile is too nar¬ 
row to explain the required gas temperature, (ii) the pres¬ 
ence of strong C IV in the FOS spectrum and (iii) the non¬ 
detection of S IV and Si HI absorption. The average chemi¬ 
cal/ionization conditions that we present in this section are in 
good agreement with those from component-by-component 
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Figure 2. The HST /WFPC2 F702W image of the Q 0122-003 field on the left, with spectroscopic redshifts of the various galaxies labelled. The top right is 
a 5 X 5" zoomed in view of the galaxy of interest, at Zg^i = 0.39853, which is classified as an She galaxy. A part of the galaxy spectrum, obtained with the 
KECK /ESI, covering Ho and N II emission lines, is shown in the lower right panel. The data are shown in black, a continuum fit is shown in red and the error 
or sigma spectrum is shown in green. 


Table 3 

Redshifts and properties of the galaxies in the Q 0122-003 field. 


Galaxy ID 

RA 

Dec 

™'F702W 

D{") 

D(kpc) 

^Kal 

i (degrees) 

$ (degrees) 

Q 0122-003G1 

01:25:29.29 

-00:05:43.13 

21.5 

14.44 

57.7±0.3 

0.25739±0.000050 

85.Ol^'^ 

3.01]'^ 

Q 0122-003G2 

01:25:28.21 

-00:05:54.64 

22.2 

9.66 

76.5±1.3 

0.95410±0.000095 

60.5li«;2 

1 T 0+45.2 
-Ld.5_;1^9 g 

Q 0122-003G3 

01:25:28.26 

-00:06:08.20 

20.4 

15.18 

78.9±0.3 

0.37874±0.000022 

68.oil's 

57.5l^'^ 

Q 0122-003G4 

01:25:27.91 

-00:05:53:17 

22.3 

14.39 

117.2±0.4 

1.07642±0.000069 

85.0l°;° 


Q 0122-003G5 

01:25:28.27 

-00:06:23.03 

19.5 

28.59 

122.2±0.8 

0.28244±0.000021 

84.3l°'i 

84.911° 

Q 0122-003G6 

01:25:27.67 

-00:05:31.39 

19.4 

30.39 

163.0±0.1 

0.39853±0.000027 

63.211;^ 

73.4lli 

Q 0122-003G7 

01:25:30.68 

-00:06:32.43 

20.3 

45.98 

239.2±0.1 

0.37923±0.000029 

56.i11°2^ 

44.7llj 


analyses. The average model parameters thus provide an ade¬ 
quate description of the high- and low-ionization gas phases. 

As illustrated in the left panel of Figure 3, the density of the 
high-ionization phase, derived using the /V(N v) to N(0 Vl) 
ratio, is —4.15 ^ lognn ^ —4.00 corresponding to an ion¬ 
ization parameter of—1.70 ^ log {7 ^ —1.85. For this given 
density, a super-solar metallicity, i.e. [X/H] > 0.3, is re¬ 
quired to reproduce the observed column densities of O VI and 
N V. This PI solution produced considerable amounts of C III 
(log N = 14.9) and C IV (log N = 15.3) but did not produce 
any other detectable low-ions (e.g. C II, Si III). 


PI models for the low-ionization gas are presented in the 
right panel of Figure 3. The density is constrained to be 
—2.75 < logriH ^ —2.40 (i.e. —3.10 ^ log U ^ — 3.45) 
from the N (Mg l) to N (Mg ll) ratio. A metallicity of 
[X/H] > —1.4 is needed to explain the observed column den¬ 
sities of Mg I and Mg II simultaneously. Such a solution does 
not produce a significant amount of N V and/or O VI. 

In Table 4 we summarize the PI model parameters. The 
large total hydrogen column densities (i.e. log A^h > 19.0) 
suggest that a significant amount of baryons is associated with 
both phases. The thicknesses (Lios) of the two phases are tens 
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Table 4 

Summary of PI model parameters. 


Phase 

log nn 
(cm-3) 

iogt; 

> [X/H| 

> log A^h 
( cm-2) 

^ -^los 
(kpc) 

High 

-4.1 

-1.75 

0.3 

19.1 

53 

Low 

-2.6 

-3.25 

-1.4 

20.2 

19 


of kpc and differ from each other by only a factor of ~ 3. 
But, the density of the low-ionization phase is ~ 30 times 
higher than that of the high-ionization phase. Most interest¬ 
ingly, the metallicity of the high-ionization phase is over an 
order of magnitude higher than the low-ionization gas phase. 
The data clearly indicate different origins of the high- and 
low-ionization gas phases. 

4. DISCUSSION 

4.1. The strong and large velocity spread O VI absorption 

We have presented a detailed analysis of an ultra-strong 
intervening O VI absorber at Zahs = 0.39853 detected in the 
HST /COS spectrum of Q 0122-003. The absorber exhibits 
logiV(Ovi) = 15.16 ± 0.04 with a kinematic spread of 
Augo = 419 km s“^. Such a strong (i.e. log A > 15.0) and 
large velocity spread (i.e. > 400 km s“^) intervening Ovi 
systems are very rare in the Milky Way halo (Wakker et al. 
2003; Savage et al. 2003), in high-velocity clouds (HVCs; 
Sembach et al. 2003), in the COM (Werk et al. 2014) and 
in the IGM (Tripp et al. 2008; Muzahid et al. 2012; Savage 
et al. 2014). For example, only 2/42 systems in the COS- 
Halos sample show log A(Ovi) > 15.0 (Tumlinson et al. 
2011b). In a blind survey of intergalactic Ovi absorbers at 
high-2:(1.9 < 2 < 3.1) Muzahid etal. (2012) have found only 
3/85 systems with Augg > 200 km s“^ with none of them 
having Avgg > 400 km s“^. The systems with large Avgo, 
in their sample, are thought to be associated with outflows 
from Lyman Break Galaxies (LBGs; see also Adelberger et al. 
2003, 2005). The only known O VI absorbers at low-z with 
Augo in excess of 400 km s“^ are the ^abs = 0.3558 sys¬ 
tem towards J 1009-1-0713 (Tumlinson et al. 2011a), the Zabs 
= 0.927 system towards PG 1206H-459 (Fox et al. 2013) and 
the Zabs = 0.227 systems towards the QSO pair Q 0107-025A 
and Q 0107-025B (Muzahid 2014). Tripp et al. (2011) have 
argued that a post-starburst outflow is responsible for the wide 
velocity spread absorption seen towards PG 1206-1-459 (but 
see Ding et al. 2003, for an alternate explanation). Tumlinson 
et al. (2011a) have suggested that Ovi arises from interfaces 
between “cool” and “hot” gas phases. Moreover, Muzahid 
(2014) has conjectured that the large velocity spread O VI ab¬ 
sorption seen through the pair of lines-of-sight is originating 
from an “ancient outflow” (see e.g. Ford et al. 2014). 

Recently, in 15 of 20 Lyman limit systems (LLS) with 
logA(Hl) > 17.3 at 2 < z < 3.5, Lehner et al. (2014) 
have reported presence of strong (e.g. average log A(0 Vl) = 
14.9 ± 0.3) and wide velocity spread i.e. 200 < Av < 400 
km s“^ O VI absorption in the ISM/CGM of high-z galaxies 
and hypothesized that these strong CGM O VI absorbers probe 
outflows of star-forming galaxies (see also Fox et al. 2007). 
We note that, a vast majority of local starburst galaxies show 
ultra-strong “down-the-barrel” O VI absorption with an aver¬ 
age column density of log A(Ovi) = 15.0 ± 0.3 (Grimes 


et al. 2009). Additionally, these O VI absorption lines are 
broad with a FHWM of - 367 ± 95 km s’^. 

In summary, intervening O VI absorbers with log A > 15.0 
and Augo > 400 km s“^, as observed in the present system, 
are extremely rare. Large-scale galaxy outflows are thought 
to be the most likely origin of such O VI absorbers whenever 
detected. The frequent occurrence of strong and large veloc¬ 
ity spread O VI absorption in local starburst galaxies further 
strengthens this idea. 

4.2. The presence of strong N V absorption 

Unlike most low-z intervening O VI absorbers, the present 
system shows an ultra-strong (log A(N v) = 14.69 ± 0.07) 
and wide-spread {Avgg = 285 km s“^) N V absorption with 
log A(Nv)/A(Ovi) = —0.47 ± 0.08. As the creation 
and destruction ionization potentials of N V are 77.5 eV and 
97.9 eV respectively, a hard radiation held at energies above 
4 Ryd is required to ionize nitrogen to its ionization 
state. Though no systematic survey exists, weak N V ab¬ 
sorption (< 10^^ cm“^) is detected in only a handful of 
O VI systems at low-z (e.g. Zabs = 0.20260 system towards 
PKS 0312-77 (Lehner et al. 2009), Zabs = 0.16716 system 
towards PKS 0405-123 (Savage et al. 2010), Zabs = 0.227 to¬ 
wards Q 0107-025A (Muzahid 2014)). In all of them O VI is 
found to be fairly strong (> lO^'^'^ cm“^). Flere, we note that 
Werk et al. (2013) have not found any N V absorption with 
log A > 14 in the COS-Halos sample. 

To our knowledge, there are only 6 other intervening sys¬ 
tems that show logA(Nv) > 14.0 and only one of them 
is at z < 1. These are the Zabs = 0.927 system towards 
PG 1206H-459 (Tripp et al. 2011), Zabs = 2.811 system to¬ 
wards Q 0528-250 (Fox et al. 2007), Zabs = 2.83437 sys¬ 
tem towards J 1343-1-5721 and Zabs = 2.47958 system to¬ 
wards Q 1603-1-3820 and Zabs = 2.18076 system towards 
Q 1217-1-499 (Lehner et al. 2014), and the Zabs = 1.5965 sys¬ 
tem towards PKS 0237-23 (Fechner & Richter 2009). Inter¬ 
estingly, all these systems show strong Ovi (log A > 14.8) 
absorption kinematically spread over > 200 km s“^ similar 
to the present system^, and all of them are LLS. 

A population of photoionized intervening N V absorbers, 
with dM/dz ^ 0.9 down to log A = 12.7, at high-z 
(1.5 < z < 2.5) is reported by Fechner & Richter (2009). 
These absorbers are predominantly weak and only one out of 
21 systems shows log A(Nv) > 14.0. As the intensity of 
the extra-galactic UV background radiation is considerably 
higher as compared to low-z due to the enhanced QSO ac¬ 
tivity at z ~ 2 (Haardt &. Madau 1996; Khaire & Srianand 
2015), detection of weak N V absorption at high-z is consis¬ 
tent with expectations for PI models. However, the presence 
of N V with log A > 14.0 is very rare at any redshift. 

Apart from the present system, the only other O VI absorber 
with log A(N v) > 14.0 at z < 1.0 is the well known system 
at Zabs = 0.927 towards PG 1206-1-459 (Churchill & Charlton 
1999; Ding et al. 2003; Tripp et al. 2011). These two sys¬ 
tems have remarkable resemblance in the absorption proper¬ 
ties of their high-ions and show log A (N v)/A(0 Vl) ^ —0.5 
(Rosenwasser et al., in preparation). Unlike Tripp et al. 
(201 1), who favored a collisional ionization scenario for N V 

^ Note that 0 VI information is not available for the z^bs = 1.5965 system 
towards PKS 0237-23. However it shows a strong CIV absorption spread 
over 300 km s“^. The strong congelation between the line spreads of CIV 
and O VI absorption found by Muzahid et al. (2012) suggests that this system 
should have a strong and wide-spread O VI absorption. 
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Figure 3. Average PI models for the high- (left) and low- (right) ionization gas. Bottom panels: The estimation of density/ionization parameters using ionic 
ratios. The horizontal dashed and dotted lines represent the observed column density ratios and their Itr uncertainties, respectively. The shaded regions indicate 
the density ranges over which model column density ratios (magenta curves) match with the observed values. The ionization parameters corresponding to 
different nu values are labelled at the top of each panel. Top panels: Metallicity required to reproduce the observed column densities of high- (O VI and N V) 
and low- (Mg I and Mg II) ions are plotted in thick curves as a function of gas density. The dotted curves correspond to Icr uncertainties in respective column 
density measurements. The JV(H l) values assumed by the models are mentioned. The horizontal dotted lines represent the average metallicities. Over an order 
of magnitude difference in metallicity between the high- and low-ionization gas phases is apparent. 


and could not estimate a high-ionization phase metallicity, 
Rosenwasser et al. have found that both O VI and N V can 
be explained under PI equilibrium with solar to super-solar 
metallicities and with a density of ^ 10“^ cm“^. As we dis¬ 
cuss next, all these properties are strikingly similar to what 
we deduce for the present system. 

Very strong N V absorption is characteristic of many intrin¬ 
sic narrow and mini-broad absorption line systems, that are 
known to be related to the AGN winds/host-galaxy because of 
a velocity within ~ 5000 km s“^ of the QSO redshift, and ev¬ 
idence of partial coverage of the QSO continuum/broad emis¬ 
sion line source (Hamann et al. 2000; Srianand et al. 2002; 
Misawa et al. 2007; Wu et al. 2010). For example, Wu et al. 
(2010) reported three 2.6 < z < 3.0 intrinsic N V absorbers 
with metallicities greater than 10 times the solar value. Strong 
O VI and C IV absorption are common in these systems as 
well. The present system is clearly not an intrinsic absorber, 
but the common properties are likely evidence of an origin in 
the high metallicity, central region of a galaxy. Both the ma¬ 
terial that feeds an AGN wind, and that which is ejected in 
a starburst outflow share a similar chemical origin in a high 
metallicity environment (Hamann & Ferland 1999). 

4.3. Origins of different gas phases 

In Section 3.3 we have demonstrated that the Ovi bear¬ 
ing gas in this absorber can be explained with PI equilibrium 
models that require a density of ^ 10“'’^'^ cm“^, a super-solar 


metallicity (i.e. [X/H] > 0.3), and a solar [N/0] ratio. The 
synthesis of N can arise from primary and secondary produc¬ 
tion via the CNO cycle. Primary production is synthesis of 
N from C produced in the core of the same star via helium 
burning and secondary production is the synthesis of N from 
C and O produced in previous generations of stars. Secondary 
N enrichment of the ISM occurs well after massive stars have 
gone Type-II supernovae and seeded the ISM with oxygen. 
For enrichment from primary synthesis, the [N/0] ratio is 
roughly —0.6 for [0/H] < —0.3. For enrichment from sec¬ 
ondary synthesis, the [N/0] ratio increases with increasing 
metallicity for [0/H] > —0.3, such that by [0/H] ~ -1-0.3, 
[N/0] ~ 0.0 i.e. the solar value (e.g. see Figure 9 of Pettini 
et al. 2008). Therefore, the assumption of a solar [N/0] ratio 
in this system is reasonable. Solar/super-solar metallicity with 
solar [N/O] ratio suggests that nitrogen in the high-ionization 
absorbing gas is predominantly produced via the secondary 
channel. Thus, the high-ionization, low density gas phase is 
likely to be originating from a region of the galaxy with a 
prolonged and high star formation rate (see also Hussain et al. 
2015). This is consistent with the host-galaxy properties as 
discussed in the next section. 

It is often claimed that strong O VI absorbers are collision- 
ally ionized (e.g. Lehner et al. 2014) as they require unrea¬ 
sonably large sizes under PI equilibrium (but see Muzahid 
2014). However, from the observed narrowness of the Ly/? 
profile, the non-detection of S IV absorption, and the presence 



















An extreme metallicity, large-scale galaxy outflow 


9 


of strong C IV absorption in the low-resolution FOS spectrum, 
we have ruled out the possibility of O VI bearing gas in this 
system being collisionally ionized (Appendix B.3 and B.4). 

For the low-ionization phase, both the overall (Sec¬ 
tion 3.3) and the component-by-component (Appendix B.2) 
approach of PI modeling suggest a metallicity of [X/H] > 
—2.0 to —1.3. The high-ionization phase has more than a 
factor of ^ 10 higher metallicity. This is in contrast to the 
Zahs = 0.927 system towards PG 1206H-459 for which Rosen- 
wasser et al. (in preparation) have found that the solar/super¬ 
solar metallicity of the low-ionization gas is indeed similar 
to the high-ionization gas. We note that in the present case 
there is no one-to-one correspondence between the high- (e.g. 
O VI) and low- (e.g. Mg ll) ionization absorption kinematics. 
But in the case of PG 1206H-459 kinematic similarity between 
Mg II and Ovi absorption is quite remarkable. Therefore, 
while for PG 1206H-459 both the high- and low-ionization gas 
phases are consistent with being part of same multiphase out¬ 
flow (e.g. Tripp et al. 2011), the low-ionization phase in the 
present system likely has a very different origin from the high- 
ionization phase. 

Using the methods of Kacprzak et al. (2010), we mea¬ 
sured the rotation curve of the host-galaxy from the Fla emis¬ 
sion line, which extends ±1" (corresponding to 3.5 kpc) 
from the galaxy center. The projected rotation velocity is 
~ 200 km s“^. We then employed the rotation model of 
Steidel et al. (2002) to examine the range of allowed veloc¬ 
ities along the line-of-sight that are consistent with extended 
galaxy rotation given the galaxy inclination, azimuthal angle, 
the quasar impact parameter, and the gas scale height, h,^. 
Following Kacprzak et al. (2010), we use h,^ = 1 Mpc (for 
a non lagging halo above the galaxy plane). The range of ve¬ 
locities along the line-of-sight that could be due to extended 
rotation from the galaxy are -1-25 < v < 4-150 km s“^. 
The components of the low-ionization gas are at u ~ 0 and 
V ~ 200 km s“^, as probed with Mg II absorption. None 
of the Mg II components is consistent with co-rotation of the 
galaxy. TTie low metallicity and kinematics of this phase per¬ 
haps suggest that it stems from recycled gas as seen in numer¬ 
ical simulation of Oppenheimer et al. (2010). Such an inter¬ 
pretation is also supported by the recent study of “bimodal” 
metallicity distribution of low-ions in low-2 LLS by Lehner 
et al. (2013). The metal-poor branch, possibly tracing cold 
accretion streams, in their sample peaks at [X/H] = —1.6 
which is consistent with what we And for this system. 

4.4. The host-galaxy 

At the epoch of our measurement, the host galaxy has a high 
SFR of 6.9 Mq yr“^ and a E* of ~ 0.4 Mq yr“^ kpc“^, the 
latter being a factor of ^ 4 above the well known threshold of 
0.1 Mq yr“^ kpc~^ for driving galactic-scale winds (Heck¬ 
man 2003). The observed galaxy orientation parameters, i.e. 
i = 63° and $ = 73°, as given in Table 3 suggest that the 
QSO line-of-sight is at ideal position to pierce one of the bi¬ 
cones of a bi-conical outflow propagating along the minor- 
axis of the galaxy (e.g. Bouche et al. 2012; Kacprzak et al. 
2012,2014; Bordoloietal. 2014b; Schroetteret al. 2015). The 
host-galaxy has a metallicity of 12-t-log(0/H) = 8.68±0.02, 
i.e. [0/H] ^ 0.0, which is roughly consistent with the metal¬ 
licity of high-ionization gas phase. All these observational 
evidences strongly indicate that the high-ionization gas, as 
probed by the ultra-strong O VI and N V absorption, is trac¬ 
ing a powerful (see the next section), active, metal-enriched, 
large-scale galactic outflow. 


Here we emphasize that the absorber is unlikely from the 
ISM of a dwarf galaxy that happens to be coincident on the 
quasar, van Zee & Haynes (2006) found that the H II regions 
of dwarf galaxies have sub-solar metallicity and [N/0] ratios 
consistent with primary N synthesis. On the other hand, van 
Zee et al. (1998) showed that the H II regions of massive spi¬ 
ral galaxies have solar to super-solar metallicity and [N/0] 
ratios consistent with secondary N synthesis. As such, we fa¬ 
vor the scenario in which the high-ionization absorbing gas 
was ejected from the host spiral galaxy we have identified at 
the absorption redshift. 

4.5. The outflow properties 

From the halo mass {Mh ~ IO^^'^M^q) and the virial radius 
(i?vir ~ 280 kpc) of the host-galaxy we estimate a circular 
velocity of Uc ~ 220 kms“^. Using the relationship between 
SFR, Vc, and the terminal outflow speed (um) by Sharma & 
Nath (2012) we estimate Vu, 230 km s“^. Here we note 
that Vc is calculated at i?200 in their models whereas we com¬ 
pute Vc at the virial radius. Nevertheless, this does not change 
our estimated Vw by more that a few percent. An outflow with 
a constant speed of 230 km s“t would require 0.7 Gyr time 
to travel the minimum projected distance of 163 kpc. This is 
on the order of the sound crossing time for the individual cool 
{T ^ 10^ K) photoionized clouds with size ~ 4 —10 kpc (see 
Appendix B.2). However, in practice the outflow could have 
been ejected with a significantly higher speed as compared 
to the terminal speed we estimate here. In that case the flow 
time could be significantly lower as compared to the sound 
crossing time so that the clouds would have enough time to 
be stable against any mechanical disturbances before they tra¬ 
verse a distance of 163 kpc. 

Under the thin-shell approximation the mass of the outflow¬ 
ing gas can be written as: Mout = 4:TrpuipCnCf NB_r^ , where 
rrip is the proton mass, p, = 1.4 accounts for the mass of he¬ 
lium, Cq (Cf) is the global (local) covering factor, and r is 
the thin-shell radius (e.g. Rupke et al. 2005). The global cov¬ 
ering factor is related to the wind’s opening angle. Here we 
assume Cq = 0.4 as observed in local starbursts (e.g. Veilleux 
et al. 2005). The local covering factor, on the other hand, is 
related to the wind’s dumpiness. For simplicity we assume 
it to be unity for the high-ionization, low density, diffuse gas. 
Using r = 163 kpc and A^h = 10^®'^ cm“^ (see Table 4) we 
estimate Mout 2 x IO^^Mq corresponding to an oxygen 
mass of Mo 10^ Mq. These are consistent with the masses 
as estimated in the CGM of ~ L* galaxies at 2 < 0.2 by Tum- 
linson et al. (201 lb). Using the wind speed of 230 km s“^, we 
further estimate the mass-flow rate of Mout ~ 54 Mq yr“^ 
and the kinetic luminosity of 77^ ~ 9 x 10^^ erg s“^ (see Ap¬ 
pendix C). These values are typical of what have been seen in 
“down-the-barrel” outflows from infrared-luminous starbursts 
at 2 < 0.5 (i.e. Rupke et al. 2005). 

To better understand how the mass outflow rate is related to 
the SFR of the host galaxy, it is customary to define the mass 
loading factor, p = Mout/SFR, which is a critical ingredient 
for numerical simulations of structure formation (e.g. Oppen¬ 
heimer et al. 2010; Dave et al. 2011b). We estimate a p of 
~ 8 for the outflow studied here. In the infrared-luminous 
starbursts sample of Rupke et al. (2005), rj is found to range 
from 0.001 to 10 but is ^ 0.1 on an average. However, we 
note that a vast majority of these starburst galaxies have SFR 
> 100 Mq yr“^ which is significantly higher as compared to 
the present galaxy with SFR of 6.9 Mq yr“^. Therefore, the 


10 


S. Muzahid et al. 


galaxy in our study is extremely efficient in entraining mass, 
in the form of outflow, for its SFR. This probably because 
of the high E*. However, it might also be possible that the 
host-galaxy had much higher SFR at the epoch of ejection 
compared to what we observe today. 

5. SUMMARY 

We have examined the physical conditions, chemical abun¬ 
dances, and energetics of a large-scale galactic outflow in the 
COM of a star-forming (SFR 6.9 Mq yr"!), sub-L* (0.5L;^) 
galaxy at z = 0.39853. Using the halo abundance matching 
technique we estimate a halo mass of Mh ~ 
and a virial radius of i?vir ^ 280 kpc for the galaxy. Along 
with Lya and Ly/3, several low- (i.e. Mg I, Mgll, Cll, N II, 
Sill, Sim) and high- (i.e. Civ, Nv, Ovi) ionization metal 
absorption lines, originating from the CGM of the galaxy, 
have been detected in the HST(COS, FOS) and ULT(UVES) 
spectra of the background quasar Q 0122-003. The QSO is lo¬ 
cated at an impact parameter of D = 163 kpc (D /i?vir = 0.6) 
with an azimuthal angle of 73° with respect to the galaxy 
major-axis. Our main findings are as follows: 

1. The low-ionization absorption components trace pho- 
toionized gas with a density of nn ^ 10“^ ® cm“^ 
(log (7 ~ —3.2) and a metallicity of [X/H] > —1.4. 
Such a low metallicity gas is consistent with being re¬ 
cycled material in the galaxy halo as predicted in nu¬ 
merical simulations. 

2. The high-ionization gas phase shows ultra-strong O VI 
with log N{0 Vl) = 15.16 ± 0.04 kinematically spread 
over Augo = 419 km s“^. Such a strong and large 
velocity spread O VI absorber is very rare at any red- 
shift. Among the known astrophysical environments 
in the local universe, only starburst galaxies show such 
strong O VI absorption (Grimes et al. 2009). 

3. Unlike most intervening Ovi systems, this absorber 
shows ultra-strong and wide velocity spread N V ab¬ 
sorption with log W(N v) = 14.69 ± 0.07 and Augg = 
285 km s“^. There is only one system at 2 < 1 (i.e. 
-Zabs = 0.927 towards PG 1206H-459, Tripp et al. 2011) 
that shows N V as strong as this. 

4. It is often claimed that strong O VI absorbers are colli- 
sionally ionized as they commonly require large sizes 
(> Mpc) under PI equilibrium (Lehner et al. 2014). 
However, for this system our PI models suggest line- 
of-sight thicknesses for individual components of ~ 
10 kpc. In fact, based on the 6-parameter constraints 
from the Ly/3 profile, the non-detection of S IV, and the 
presence of strong C IV (Wr = 1.7 A) in the FOS spec¬ 
trum, we have ruled out CIE and non-CIE models. 

5. The high-ions, O VI and N V, can be well explained as 
arising in a low-density (nn ~ 10“^'^ cm“^, log U ~ 
— 1.6) photoionized gas with super-solar metallicity 
([X/H] > 0.3) and with solar [N/0] ratio. The super¬ 
solar metallicity with a solar [N/O] ratio implies that 
nitrogen in this system is predominantly produced via 
secondary synthesis. Thus the high-ionization, low- 
density gas is presumably stemming from regions of the 
host-galaxy that have sustained a high star-formation 
rate for a prolonged period. 


6. The measured E* of ~ 0.4 Mq yr“^ kpc“^ for the host 
galaxy is a factor of ~ 4 higher than the threshold of 
^0.1 Mq yr“^ kpc“^ for driving galactic-scale wind 
as seen in local starbursts. The azimuthal angle of 73° 
suggests that the QSO line-of-sight is in ideal position 
to pierce the bi-conical outflow propagating along the 
minor-axis of the host-galaxy. 

7. Assuming the high-ionization gas is in outflow, we es¬ 
timate a outflow mass of Mout 2 x 10^° Mq, mass 
outflow rate of Mout 54 MQyr“^, and kinetic lu¬ 
minosity of i/fc ^ 9 X 10"^^ erg s“^ using a thin-shell 
model. These values are consistent with that of “down- 
the-barrel” outflows from infrared-luminous starbursts 
at z < 0.5 (Rupke et al. 2005). 

8. We estimate a mass loading factor of p ^ 8. This is 
among the highest values as seen in infrared-luminous 
starbursts at low-z. However, we note that the SER of 
the host-galaxy is significantly lower (by a factor of ten 
or more) than the starburst galaxies (Rupke et al. 2005). 
This indicates that the host-galaxy is highly efficient in 
entraining mass onto the CGM in the form of outflow. 

Einally, we emphasize that such powerful, large-scale, 
metal-rich outflows are the essential component by means of 
which a galaxy imparts sufficient mechanical and chemical 
feedbacks that regulates star formation and hence evolution 
of a galaxy. Einding and analyzing more such unique absorp¬ 
tion systems is crucial for a comprehensive understanding of 
galaxy feedback mechanisms and reinforce galaxy evolution 
theory with useful observational constraints. 
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APPENDIX 

A. ESTIMATION OE N{H l) IN DIEEERENT LOW- 
AND HIGH-IONIZATION COMPONENTS 

In this section we describe how A(Hl) values in differ¬ 
ent low- and high-ionization absorption components are es¬ 
timated. Determining N (H l) component-by-component in a 
multiphase, multicomponent absorption system is extremely 
challenging. The presence of unsaturated higher-order Lyman 
series lines is very helpful in order to disentangle the contribu¬ 
tions from high- and low-ionization phases (see e.g. Muzahid 
2014). Our present COS spectrum only covers Lya and Ly/3 
absorption for this absorber. We, therefore, try to estimate 
conservative upper limits on A(H l) in different low- (L-1 to 
L-3) and high- (H-1 to H-5) ionization components. 
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Figure Al. Estimation of upper limits on A^(Hl) in different low- (left panel) and high- (middle panel) ionization absorption components. Individual profiles 
are shown in (green) dotted curves. The resultant profiles are shown in solid curves. The error in each pixel is shown as a grey error-bar. The vertical dotted 
lines in the left panel are line centroids of the Mg II components. The small ticks in the middle panel are line centroids of the O VI components. Both low- and 
high-ionization components ai‘e shown in the right panel. Contributions from low- and high-ionization components are shown in (blue) dotted and (red) dashed 
curves respectively. The resultant solid (magenta) curve explains both the Lya and Ly/3 profiles satisfactorily. 


To estimate the maximum allowed N{H l) associated with 
the low-ionization components, we use the Mg II compo¬ 
nent structure and assume pure non-thermal line broadening, 
i.e. 6(Hl) = 6(Mgll)"^. This gives the lowest possible b- 
parameters for H I in each component. We then use VPFIT to 
generate synthetic profiles corresponding to different V(H l) 
values. We start with an initial value of log TV(Hi) = 14.0 
and gradually increase the V(Hl) value until the synthetic 
profile just exceeds the observed Lya and/or Ly/3 profiles. 
The use of the lowest possible 6(H l) values ensures that the 
N{H l) values we report in Table 1 are indeed conservative. 
The left panel of Figure Al shows the synthetic profiles cor¬ 
responding to the maximum allowed 7V(H l) values we esti¬ 
mate for the low-ionization components. It is apparent from 
the bottom panel that the wings of the Lya absorption are re¬ 
produced adequately. Note that because of the assumed low 
6-parameters, V(Hl) value as high as ^ 10^® ® cm“^ does 
not produce any significant Ly/3 optical depth. Therefore, al¬ 
most all the Ly^ absorption and the remaining absorption seen 
in Lya must arise from the high-ionization components. As 
mentioned in the footnote of Table 1 our V(H l) limits remain 
unaltered for 6-values in the range 14 - 20 km s“^ which is 
expected from the PI model temperatures. 

To estimate the maximum allowed 7V(H l) corresponding 
to the high-ionization components, we use the O VI compo¬ 
nent structure and again assume pure non-thermal line broad¬ 
ening, i.e. 6(Hl) = 6(0 Vl). Similar to the low-ionization 
components, the assumption of pure non-thermal broadening 
ensures that the V(Hl) upper limits for the high-ionization 
components are conservative as well. We fit the remaining 
absorption (i.e. after accounting for the contribution from 
the low-ionization components) using VPFIT by allowing the 
code to vary V(Hl). The maximum allowed V(Hl) for the 
high-ionization components are given in Table 2. In the mid¬ 
dle panel of Figure Al we show the corresponding profiles as 
dotted curves. The solid curve represents the resultant profile. 
We emphasize here that due to the presence of low-ionization 
components at similar velocities, the 7V(Hl) limits for the 
components H-2 and H-4 are uncertain. However, the N{H l) 
limits in the rest of the high-ionization components are robust. 
For the H-1 and H-5 components, the N (H l) values are con¬ 
strained adequately by the blue and red wings of the Ly/3 pro- 

At(H l) values can be significantly lower if one assumes pure thermal 
broadening, i.e. 6(Hl) = \/24 fe(Mgll). However, the wings of the Lya 
absorption cannot be explained in such a scenario. 


file, respectively. Moreover, the non-black absorption seen in 
the Ly/3 at u ^ -1-100 km s“^ provides a robust constraint on 
the iV(H l) for the H-3 component. A free fit to this part of 
the Ly/3 absorption, without considering any contributions it 
might have from the adjacent absorption, leads to a solution 
with6(Hl) = 28±4kms-i and log TV (Hi) = 15.10 ±0.10. 
This free fit 6-parameter is consistent with our assumption that 
6(Hl) = 6(0 Vl) within la allowed uncertainty. Note that 
no low-ions are detected in any of these components (i.e. in 
H-1, H-3, and H-5). Therefore, the observed Hi absorption 
is entirely due to the high-ionization phase. Next, we assume 
that the Ly/3 profile is completely due to the high-ionization 
absorption components with 6(H l) = 6(0 Vl). Under such an 
assumption the column densities we obtain are very similar to 
the upper limits we derived in the presence of low-ionization 
components (see Table 2). In the right panel of Figure Al we 
show the resultant (high- and low-ionization) profiles which 
explain both the Lya and Ly/3 profiles adequately. 

B. IONIZATION MODELS: 

COMPONENT-BY-COMPONENT 

B. 1. PI models for the high-ions 

Here, we model the three high-ionization components 
(H-3, H-4, and H-5) in which both O VI and N V are posi¬ 
tively detected and their column densities are well measured. 
As demonstrated in the previous section, the /V(Hl) limits 
in H-3 and H-5 are also robust. Eor each of these compo¬ 
nents we compute ionization fractions for different species 
as a function of density for a given input metallicity using 
the corresponding N (H l) limits as the stopping criterion for 
CLOUDY. In order to explore possible dependence of ioniza¬ 
tion fractions on the input metallicity, models were run for 
a series of different input metallicities (e.g. [X/H] = —1.0, 
0.0, 0.3, 0.7, and 1.0). We, however, find that the ionization 
fractions do not show any significant dependence on the input 
metallicity in the optically thin regime for [X/H] < 0.3 for 
the entire density range of interest. Eor the components with 
derived metallicity [X/H] > 0.3, we use the CLOUDY model 
for which input metallicity is close to the derived metallic¬ 
ity. Using the model predicted ionization fractions and the 
observed column densities we determine loci of O VI and N V 
in the density-metallicity (lognn — [X/H]) plane. 

In different panels of Figure A2, these loci of O VI and N V 
are shown. The curve corresponding to each ion (O VI or N v) 
represents the required metallicity to reproduce the observed 
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Figure A2. PI models for the high-ionization absoiption components, H-3 (left panel), (middle panel), and H-5 (right panel). Metallicity required to 
reproduce the observed column densities of individual high-ions (O VI and N V) are plotted as a function of gas density. The Al(Hl) values assumed by the 
models ai'e mentioned in the plot (see Table 2). The ionization parameters (log U) coiresponding to different nji values ai'e labelled at the top of each panel. The 
shaded regions represent the range in log njj (or log U) over which N{0 Vl) and Af(N v) match with our measurements. The horizontal dotted line represents 
the lower limit in metallicity indicated by the arrow. 




log U = -2.5 -3 -3.5 -4 



Figure A3. Same as Figure A2 but for the low-ionization absorption components, L-1 (left panel), L-2 (middle panel), and L-3 (right panel). The measured 
Af(Mg II) and Af(Mg l) values along with the N (H l) limits are used to constrain the model parameters. As Mg I is not detected in L-1, N{C It) is used instead. 


column density as a function of log riH (or log U). The range 
in [X/H] and in lognn for which both the curves overlap, 
represent a PI equilibrium solution since both 7V(Ovi) and 
iV(N v) are simultaneously reproduced. 

In Table B1 we summarize the PI model solutions for the 
three high-ionization components. Densities as estimated for 
different components are consistent within ^0.4 dex with a 
median value of nn ~ 10“^'^ cm~^ (i.e. \ogU = —1.65). 
The most interesting aspect of our models is to note the con¬ 
servative lower limits on gas phase metallicities. All three 
components show super-solar metallicities with [X/H] > 0.3. 
For the H-3 and H-5 components, in which the N (H l) limits 
are well constrained, the lower limits turn out to be [X/H] > 
0.6 and [X/H] > 0.7, respectively. These metallicities are 
very high and very rare for any intervening absorption sys¬ 
tem. Such super-solar metallicities, however, is common in 
intrinsic QSO absorbers (e.g. Petitjean et al. 1994; Wu et al. 
2010; Muzahid et al. 2013) that stem from the black hole ac¬ 
cretion disk-wind. The total hydrogen column density is in 
the range log A^h = 17.8 — 18.6 and the line-of-sight thick¬ 
ness range from 4-10 kpc. Note that the PI solutions for 
the high-ionization components did not produce any consid¬ 
erable column densities for the low-ions (e.g. CII, Si II, Si III) 
but produced strong C III, C IV, and moderate Si IV column 
densities. Finally, we find that the PI modeling for the H-1 
and H-2 components, in which N V is not detected, suggest 
logriH < —4.1 and [X/H] < 0.0. 


B.2. PI models for the low-ions 

In Section 3.1 we have demonstrated that the kinematics 
of O VI and Mg II absorption are significantly different which 
indicate the multiphase structure of the absorber. In the previ¬ 
ous section, we have shown that the gas phase that gives rise to 
O VI and N V under PI equilibrium cannot produce any con¬ 
siderable low-ions column densities. This further reinforces 
the idea of multiphase nature of the absorbing gas. To con¬ 
strain the ionization parameters and chemical abundances in 
the low-ionization absorption components we have used the 
observed column densities of A^(Mg l) and A^(Mg ll) and the 
upper limits on A^(Hl). Note that, as Mg I is not detected, 
N{C ll) is used for constraining the model parameters for the 
L-1 component. 

PI models for the low-ionization components L-1, L-2, and 
L-3 are shown in Figure A3 and the model parameters are 
summarized in Table B2. The densities in all these compo¬ 
nents are very similar (i.e. nu ~ j^q-2.5±o.i The 

metallicity in the L-1 component is [X/H] > —2, but it is 
roughly ~ 0.7 dex higher for the L-2 and L-3 components. It 
is interesting to note that while for the high-ionization com¬ 
ponents with detected N V, data do not allow a metallicity of 
[X/H] < 0.3, the low-ionization components can have signif¬ 
icantly lower metallicities. However the metallicity need not 
necessarily be this low. Here, we have placed a very conser¬ 
vative lower limit on metallicity. The total hydrogen column 
densities in these components (i.e. log A^h ~ 19-7 — 19.8) 
are > 10 times higher than what we have found for the 
high-ionization components (see Table Bl). The line-of-sight 
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Table B1 

Photoionization model parameters for the high-ionization absorption components. 


Component lognj^ log U > [X/H]“ > logA^n > Lios logT Predicted log N (cm ^) 

('Vnr'l ('K''! N TT Si TT Si TTT P TT P TTT 



(cm-'^) 



(cm- 2 ) 

(kpc) 

(K) 

N II 

Sin 

Sim 

Cll 

Cm 

Civ 

Siv 

( 1 ) 

( 2 ) 

(3) 

3) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

TTOT 

( 11 ) 

Tnr 

(13) 

iT?r 

H-3 

1-4.25, -4.15), (-4.20) 

1-1.70, -1.60] 

0.60, ( 0 . 8 ) 

18.1 

6 

3.7 

11.8 

10.8 

12.2 

12.9 

14.7 

14.8 

13.2 

H^ 

[-4.05, -3.85], (-3.95) 

[- 2 . 00 , -I. 80 ] 

0.30, (0.6) 

18.6 

10 

4.0 

12.3 

11.4 

12.8 

13.3 

15.0 

15.1 

13.7 

H-5 

[-4.35, - 4 . 25 ], (-4.30) 

[-1.60, -1.50] 

0.70, (0.9) 

17.8 

4 

4.5 

11.6 

10.9 

12.1 

12.9 

14.6 

14.6 

13.1 


Note- “Values in the parenthesis are used to compute column densities in columns 8 through 14. 


Table B2 

Photoionization model parameters for the low-ionization absorption components. 


Component log njj log U > [X/H]“ > log A^jj > Li^s log T Predicted log N (cm ^ 



(cm-'^) 



(cm-2) 

(kpc) 

(K) 

N II 

Sill 

Sim 

Cll 

Cm 

Call 

Fell 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

L-1 

-2.40 

-3.45 

-2.00 


4 

4.1 

13.3 

12.9 

12.8 

13.9“ 

13.8 

11.1 

12.6 

L-2 

-2.60 

-3.25 

-1.30 

19.8 

8 

4.1 

14.1 

13.6 

13.8 

14.5 

14.8 

11.7 

13.1 

L-3 

-2.50 

-3.35 

-1.40 

19.7 

5 

4.1 

13.9 

13.4 

13.5 

14.3 

14.5 

11.6 

13.0 


Note- “C II is used for constraining the model. 


thicknesses of 4-8 kpc are, however, comparable to those of 
the high-ionization components. Note that our low-ionization 
phase solutions do not produce any appreciable amount of 
Si IV, C IV or any other higher ions (i.e. logiV < 13.0). The 
model predicted column densities for different low-ions are 
presented in columns 8 through 14 of Table B2. We use these 
column densities to generate synthetic profiles using VPFIT 
assuming pure non-thermal broadening. The synthetic pro¬ 
files are shown in magenta curves in Figure 1. Note that 
these synthetic profiles explain the observed data adequately 
for C II, N II, Si II, and Si III. Our models, however, predict 
0.3 and 0.5 dex higher Ca II column densities in L-2 and L-4 
components, respectively. It is known that Ca easily depletes 
on to dust. To be consistent with the data our PI model so¬ 
lutions would require a maximum of 0.5 dex depletion for 
calcium. Further, we note that our model predicted N{¥q ll) 
values are 0.4 dex higher as compared to the measured val¬ 
ues in the L-2 and L-3 components. This can be understood 
in terms of a-enhancement. Note that, although we have 
built very reasonable PI models, the metallicity of the low- 
ionization components (L-2 and L-3 in particular) are not 
very well constrained due to the lack of higher order Lyman 
series lines. Spectral coverage of higher order Lyman-series 
lines down to the Lyman limit will be very useful in order to 
estimate metallicities in different low-ionization components 
accurately. Nevertheless, the metallicity of the L-1 compo¬ 
nent cannot be considerably higher than [X/H] = —2.0, in 
order to explain the red wing of the Lya profile. 

B.3. CIE models for the high-ions 

In this section we will explore the viability of CIE models 
for the high-ions. CIE models for the H-3, H^, and H-5 
components, in which N V is positively detected, are shown 
in Eigure B1 . We adopt the CIE models of Gnat & Sternberg 
(2007). In their models, equilibrium and non-equilibrium 
cooling efficiencies and ionization states for low-density ra- 
diatively cooling gas are computed under the assumption that 
the gas is dust-free, optically thin, and subject to no external 
radiation field. In each panels of Eigure B1 , we plot the metal¬ 
licity required to produce the observed N (O Vl) and A^(N v) 
as a function of gas temperature (log T) under CIE conditions. 


It is apparent from the figure that the observed N{0 Vl) and 
A^(N v) are simultaneously explained for gas temperature in 
the range, logT > 5.3 — 5.4. The lower limit on metallicity 
for the H-3 and H-5 components, in which A^(H l) values are 
robustly constrained, turns out to be [X/H] > —1.2 and for 
the H-4 component, [X/H] > —1.7. Therefore, in contrast to 
the PI models, the gas-phase metallicity can be significantly 
lower under CIE conditions provided the gas temperature is 
^ 5.3. Such a high temperature corresponds to a Hi 
Doppler parameter of fc(Hl) > 57 km s“^. However, this 
large 6(H l) value, as suggested by the CIE models, is not sup¬ 
ported by the COS data. As mentioned earlier in Appendix A 
a free fit to the non-black absorption seen in the Ly/3 pro¬ 
file at ti ^ -1-100 km s“^ (i.e. the H-3 component) gives 
6(Hl) = 28 ± 4 km s“^. This is half of the 6(Hl) value re¬ 
quired by the CIE models. Eurthermore, we find that the CIE 
solutions do not produce any considerable W(C iv) in any of 
the components. The maximum A^(Clv) predicted by our 
CIE model is log N = 13.9. However, in the low resolution 
HST /EOS spectrum, Chen et al. (2001) have detected a very 
strong C IV absorption from this system with 14/. = 1.70 A^. 
Eor comparison, the measured rest-frame equivalent width of 
O VI A1031 absorption is Wr = 1.02 ± 0.02 A. Therefore, the 
C IV is even stronger than the O VI in this system, since the / A 
value for the C IV A1548 transition is > 2 times higher than 
that of Ovi A1031. Nonetheless, CIE models fail to repro¬ 
duce strong enough C IV along with N V and O VI. These facts 
indicate that a pure CIE model is not viable. Photoionization 
models (Appendix B.l), however, can reproduce strong C IV 
along with N V and O VI from a single gas phase. 

B.4. Non-CIE models for the high ions 

In this section, we investigate if non-CIE models can ex¬ 
plain the high-ionization absorption components. Note that 
non-equilibrium effects become important in metal-rich envi¬ 
ronments when cooling efficiencies are increased and gas can 
cool faster than it recombines. Here, we consider two dif- 

^ Note that we attempted to observed the CIV absorption using the G185M 
grating. However, the resultant spectrum fails to put any useful constraints to 
our models due to poor S/N (< 2 per resolution element). 
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Figure Bl. CIE models for the high-ionization absorption components, H-3 (left panel), H-4 (middle panel), and H-5 (right panel). In each panel, different 
curves represent required metallicities to reproduce observed column densities for different high-ions (O VI and N V) as a function of gas temperature. The 
A^(H l) values assumed by the models are mentioned in the plot (see Table 2). The region in which the two curves overlap represents a CIE solution. It is evident 
that for all three components the required gas temperature is log T > 5.3, in order to reproduce the observed A^(0 Vl) and A^(N v) simultaneously. The an'ow 
indicates the lower limit on metallicity in each panel. 





Figure B2. Non-CIE isobaric cooling models for the high-ionization absorption components. Everything else is same as in Figure Bl. It is evident that the 
observed A^(0 Vl) and A^(N v) can be reproduced simultaneously for gas temperature in the range log T ~ 5.1 — 5.2. 





Figure B3. Non-CIE isochoric cooling models for the high-ionization absorption components. Everything else is same as in Figure Bl. It is apparent that the 
observed N{0 Vl) and A^(N v) are not reproduced simultaneously for any of the components for the entire parameter space indicating no non-CIE solution. 


ferent non-CIE models of Gnat & Sternberg (2007), namely, 
isochoric (constant density) and isobaric (constant pressure) 
models. As discussed by the authors, isobaric cooling is 
less rapid compared to isochoric cooling and departures from 
equilibrium are somewhat smaller. Moreover, cooling that is 
initially isobaric eventually transforms to isochoric cooling 
(see their Figure 1). 

In Figure B2 we show non-CIE isobaric cooling models for 
the H-3, H-4, and H-5 components. Clearly, these models 
suggest a gas temperature (e.g. logT ~ 5.1 — 5.2) which 
is only slightly lower than that predicted by the pure CIE 
models (see Figure Bl). Note that for the H-3 component, 
the observed 7V(0 Vl)/V(N v) ratio is consistent for a large 
range in gas temperature (i.e. logT ^ 4.7 — 5.3). The 
lower limit on metallicity for logT ~ 5.3 is [X/H] > —1.0. 


The required metallicity sharply increases as the gas temper¬ 
ature decreases. For example, any temperature lower than 
10^'^ K would require a super-solar metallicity. We note that 
for logT < 5.1 our model overproduces N{S iv) by ~ 0.4 
dex. Moreover, a temperature of logT < 4.8 overproduces 
7V(Si III) by ^ 2 dex. We, thus, rule out the possibility of 
a gas temperature lower than 10®'^ K. Unlike CIE models, 
an isobaric cooling model solution with [X/H] ~ 0.0 and 
logT ^ 5.1 can produce a large CiV column density (e.g. 
log N > 15.1). Nonetheless, as discussed in the previous sec¬ 
tion, the 5-parameter required for such a high temperature (i.e. 
5(H l) > 46 km s“^) is not supported by the Fy/3 profile. 

In Figure B3 we show non-CIE isochoric cooling models 
for the H-3, H—4, and H-5 components. Note that we do 
not find any possible solution that could explain the observed 
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A^(Ovi) and A^(Nv) simultaneously within a la measure¬ 
ment uncertainty. This clearly indicates that a non-CIE iso- 
choric cooling model is not feasible. However, within 2a 
allowed uncertainties we can have solutions which can be, 
again, ruled out by the arguments that were used previously 
to rule out the isobaric models. 


C. EQUATIONS EOR OUTELOW ENERGETICS 

In all the following equations: r is the thin-shell radius, Ar 
is the cloud thickness, is the outflow speed, T is the cloud 
temperature, A^h is the total hydrogen column density, Co 
is the global covering factor, Cf is the local covering factor, 
fflow is the flow time, Cound is the sound crossing time, Mout 
is the mass carried by the thin-shell, Mont is the mass-flow 
rate, and Ek is the kinetic luminosity. All other symbols have 
their usual meanings. We refer the reader Rupke et al. (2005) 
for detailed derivations of the Equations C3, C4, and C5. 


iflow = 5 X 10 * yrs 


too kpcy V200kms“l 


tsound — 8 X 10 yrs 


Ar 


10 kpc / \ 10^ K 


(Cl) 

(C2) 


Mout = 4:TTfj,mpCnCfNfs^r^ = 5.6 x 10 ^°Mq 


Cn„ 


Nil 


1020 ( 


100 kpc 


(C3) 


Mout = SirfirripCnCfNiirVia = 230 Mq yr ^ 
Nn 


0.4 ^ 


in20 , 




N (C4) 


Ek = 


X 



iTTkimpC^CfNfirVw^ = 2.9 X 10"*^ ergs ^ 


( ) 

( ^ 

{ V-tn \ 

V 10^*^ cm“2 J 

VlOO kpcy 

V200 kms-12 
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